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Because of the peculiar coupling of spatial and spin degrees of freedom effected by the Rashba 
spin-orbit interaction, lateral confinement of a two dimensional electronic system leads to a finite 
transverse spin polarization near the longitudinal edges of a current carrying quantum wire. The 
sign of this component of the polarization is opposite at the two edges and can be reversed upon 
inversion of the current. Interestingly for small spin orbit coupling this is the largest contribution 
to the total polarization, its magnitude being of second order in the coupling constant. As a 
consequence this phenomenon cannot be revealed in lowest order perturbative approaches. An in 
plane spin polarization component is also present that is perpendicular to the current. Within the 
same model this component would be also present in the bulk. On the other hand while in the latter 
case its magnitude is linear in the coupling constant, we find that it only represents a third order 
effect in the wire geometry. Our results are consistent with a general rigorous argument on the 
parity of the components of the spin polarization with respect to the sign of the spin orbit coupling 
constant. 

PACS numbers: 72.25.Dc, 85.75.-d, 71.70Ej, 73.63.Nm 



-a 
c 

o 
o 



> 
in 



o 

o 

-i— > 



i 

C 

o 
o 



X 



We consider in this paper the effects of the lateral 
confinement on the properties of a two-dimensional (2D) 
electron liquid in the presence of linear Rashba spin or- 
bit. The basic idea underlying our work stems from the 
simple observation that in order to construct a localized 
wavefunction one must superimpose plane waves with op- 
posite components of the wave vector along the direction 
of localization. Because spin-orbit forces the latter to 
have opposite spin components along the unrestricted di- 
rection, it is clear that in the direction of confinement the 
electronic states must display a correspondingly inhomo- 
geneous spin polarization whose direction also changes 
in space. In the particular case in which a 2D electronic 
paramagnetic system is laterally confined to produce a 
quantum wire, the superposition involves spins that are 
parallel to the plane of motion and therefore combine to 
give an out of plane polarization component. Since the 
sign of this component of the magnetization is opposite 
for states with opposite momenta along the wire, it is 
clear that a net polarization can be established when the 
occupation of such states is unbalanced as in the case in 
which a net steady state current is forced through the sys- 
tem. As we will show this net perpendicular static spin 
polarization is most pronounced near the lateral bound- 
aries and displays an oscillatory behavior. Inverting the 
direction of the current simply flips the polarization. No 
spin current is involved in the problem. 

In order to demonstrate and exemplify this idea, we 
present here a complete analysis of the simplest possi- 
ble model hamiltonian describing a system of non inter- 
acting electrons subject to a homogeneous neutralizing 
background, a linear Rashba spin orbit coupling and a 
lateral confining potential^ The effect of impurities is 
neglected although impurity scattering in the channel 
will not change our main conclusions. Although more 
realistic model potentials could be readily employed we 



have moreover chosen to describe the latter by means 
of infinite potential barriers, an assumption that allows 
us to present a transparent and explicit discussion while 
at the same time affording us a full description of the 
physical phenomenon. A great advantage of our model 
is that an analytic perturbative expansion in the spin- 
orbit coupling constant can be developed, a procedure 
that shows how the effect would not be revealed by low- 
est order perturbation theories. A generalization to more 
general smooth gating potentials is readily implemented. 

That a steady state current can lead to a spin polar- 
ization in a 2D electron liquid subject to Rashba spin 
orbit coupling is hardly unexpected. It is in fact well 
known that a current forced along the x direction of a 
bulk 2D electronic system will be accompanied by an in 
plane polarization along y£&&£& For small spin orbit 
this polarization is linear in the coupling constant. 

As we shall see however, while this component of the 
polarization is still present, confinement leads to a spin 
polarization that, for small spin orbit, is mostly out of 
plane. Interestingly, as already stated, this polarization 
can be flipped simply by inverting the current. 

The magnitude of the effect strongly depends on the 
strength of the spin orbit coupling. A quick order of 
magnitude estimate can be obtained by inspecting the 
caption of our Figure [7] while making reference to our 
model hamiltonian (Eq. (TJ)). 

While the problem can be completely solved by numer- 
ical means we will show that a careful study of the per- 
turbation expansion in the spin orbit coupling constant 
is crucial in understanding the interplay of the various 
components of the polarization. 

The paper is organized as follows: Section[T]defines our 
model and details our formulation; Section |TT] describes 
the numerical and perturbative solutions to the single 
particle problem providing a discussion of various prop- 
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erties of the eigenfunctions and eigenenergies. Section UTT1 
is devoted to the study of the total spin polarization in 
the presence of a steady state current. Many of our most 
relevant results are presented in this Section. Finally 
Section HVl provides a discussion of the most interesting 
physics including that of the crossover of the in plane 
component of the polarization and our conclusions. 



FORMULATION OF THE PROBLEM AND 
BOUNDARY CONDITIONS 





The model Hamiltonian we consider is the following: 

A2 



H = 2m + a ^ x P y ~ + V ^ 



(1) 



where the motion is restricted to the x,y plane and the 
lateral confining potential is given by: 



V(y) = 



o 

DC 



for \y\ < W/2 
for \y\ > W/2 



(2) 



Since the system is translationally invariant along the 
x direction, the eigenfunctions of the problem can be 
constructed by linear superposition of free space (V(y) = 
0) eigenstates with definite values of k x . Furthermore 
the reflection symmetry x — > — x allows one to generally 
assume k x > 0. 

For reference, we recall here that the free space electron 
eigenfunctions and eigenenergies are given by 



<Pk,±( x >v) 



i(k x x+k y y) 



±1 



*2) 



2m 



=F ah J k 2 . 



U2 

K y 



(3) 



(4) 



where L is the linear size of the system and 0k is the angle 
between the direction of the wave vector and the x-axis. 
Interestingly the Rashba spin-orbit coupling forces each 
plane wave state to have its own distinct orientation of 
the spin quantization axis. This direction lies in the x — y 
plane and makes an angle of ^ with k. These states form 
two split bands which are characterized by opposite chi- 
rality in the sense that states with the same wave vector 
have opposite spin directions in the two bands. 7,8 More- 
over, because of time reversal, in the presence of Rashba 
spin orbit, states within the same band with opposite 
wave vectors have opposite spins. The corresponding sur- 
faces (lines) at constant energy e are concentric circles 
with radii K±. These are given as the positive solutions 
of the equation 



h 2 K 2 ± 

2m 



=F H a K± . 



(5) 



Constant energy lines are schematically displayed in Fig- 
ure [1] 



FIG. 1: Possible free space states that can be superimposed 
to obtain confined eigenstates. On right we show a situation 
for which the two states of the upper branch have Reffcy] = 0. 



At a given energy e, there are only four free space eigen- 
states that can be combined to form solutions satisfying 
the zero boundary condition of the inhomogeneous prob- 
lem. This is depicted in Figure [TJ for each spin subband 
± we have two states with opposite value of k y . The four 
possible values of k y are ±fc^ , with 



fc* — \l K± — k x i 



(6) 



where K± are given in Eq. (O above. 

Since we are interested in spatially undamped solu- 
tions, we always take k x < K + . However, it will nec- 
essary to include the case K- < k x < K +1 where 
ky = i\ky\ is purely imaginary. We also define the angles 
4>± by means of: 



arctan ■ 



If k y — i\k y \ we have that also 0_ 
and from the above formula we get 



(7) 

is purely imaginary 

\k~ I 

d>- = i arctanh -r^— . 



It is convenient to consider the following two couples of 
(unnormalized) states, which are also eigenstates of the 
reflection symmetry y — > —y: 



ip± T (e,k x ,y) = 



cos(fc±y- ^) 



±i cos (kyy 

sin (tfy- ^) 
±i sin (k^y + ^) 




(8) 
(9) 



The first ± subscript refers to the chirality of the spin 
subband while the second ± subscript refers to the reflec- 
tion operation y — > —y. The eigenstates in the confined 
system are then expressed in the following way: 

<p+(k x ,y) = c++tP++ + c"+^_+ , (10) 
<p~(k x> y) = C+-V+- +c""^__ , (11) 

where the superscript in <p„ refers to the reflection parity, 
and not to the chirality. As one would expect, the solu- 
tions of the confined system are an admixture of states 



FIG. 2: Energy spectrum of a W = 0.1 /j,m wide quantum 
wire. For clarity the exact energy bands are only shown (dots) 
for discrete values of kx, while the perturbative results from 
Eq. ((23} are also plotted (solid lines) for comparison. The 
parameters values are as follows: ha — 2 10 -9 meVm and 
m — 0.3mo. The horizontal dashed line marks a Fermi level 
of 6 meV and is significant in relation to Figures [5] IS] and [7] 
The corresponding electron density is 0.67 10 12 cm -2 . 




FIG. 3: Same as in Figure [5] but with a larger value of the 
spin orbit coupling: ha — 5 10 -9 meVm. Strong deviations 
from perturbation theory are due to anti-crossing of bands. 



II. SOLUTION OF THE PROBLEM 



Exact solution at k x = 



with different chirality. A spin subband index will be 
introduced as an approximate quantum number in the 
case of small spin orbit in a following Section where we 
present a perturbative treatment of the problem. 

Imposing the boundary condition at y — W/2 for the 
(fin leads to the following conditions: 



' CO s(5^±) sin(i^) 
UCO s(^±) ^sin(S^) 



'sin(^) . , , . 
. . ,k+W+<j>+. . ,k-w+$-, 
jsin(- — g J ~ «cos(— — g ) 



,KW- 
cos ( — 




At fixed k x , the the determinants of (jT2J) and (fT3")) are 
functions of e, through the implicit dependence of k^ 
and (j)±. The zeros of these determinants provide two 
sets of discrete energies e^(k x ), corresponding to states 
of opposite parity. 

If one is not interested in separating states with dif- 
ferent reflection parity, one can combine the conditions 
ensuing from (|12p and (|13[) into the compact equation: 



( 1 — cos ky W cos k y W) sin <f) + sin </>_ 

+ (1 — cos0 + cos</>_) sinkyWsmkyW = , 



(14) 



which gives the whole spectrum for a given k x . The eigen- 
functions, dropping the reflection parity index, are hence- 
forth quite generally denoted as tp u (k x ,y). 

We finally notice that for values of e and k x such that 
ky is purely imaginary, the structure of the above formu- 
las still holds, and one need only substitute hyperbolic for 
trigonometric functions. 



A case in which the solutions can be expressed in sim- 
ple closed form is that of k x =0. For this particular value 
of k x it is immediate to see that Eq. ([T]) admits solutions 
which are eigenstates of a x . Assuming the following form 
for the wavefunctions 



^{y) is found to satisfy: 



2m dy 2 



V(y) 



l 



<%) = e<%) 



(15) 



(16) 



This result is quite general holding for an arbitrary form 
of the confining potential V(y). For the case of an infinite 
well we obtain: 



h 2 /7m\ : 
2m~ \w) 



1 



-ma 



(17) 



where n is a positive integer. For every eigenvalue there 
is an additional twofold spin degeneracy. 



B. Numerical results 

We reproduce in Figures numerical results, calcu- 
lated with the methods of the previous Section. In Figure 
[5] and [3] we present two examples of the energy spectrum 
as function of k x for different values of the spin-orbit 
coupling a. 

The properties of the eigenfunctions are perhaps more 
interesting and can be best exemplified by the corre- 
sponding number and spin polarization densities. Figures 
SI [5] and [S] refer to these properties of the eigenfunctions, 
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and have been obtained with the same parameters of Fig- 
ure [5] and an occupation characterized by a Fermi energy 
ep = 6 meV. In particular, we show there for each state 
the relative number density N v {k x ,y) and spin polariza- 
tion density Vuikx^y). These quantities are defined as 
follows: 

N v (k x ,y) = (ip v (k x ,y')\5(y -y')\ip v (k x ,y')) , (18) 



V v {k x ,y) = (tp v {k x ,y')\aS(y -y')\ip v (k x ,y')) . (19) 

From simple symmetry considerations, it follows that 
the polarization is vanishing in the x direction, along the 
wire, and states with the same energy and opposite values 
of k x have the same density but opposite spin polariza- 
tion. As it will be discussed in a following Section, this 
fact allows in principle to produce a net spin polariza- 
tion along the y and z direction, by driving an electrical 
current along the wire. 



C. Perturbation theory 

Although the numerical analysis provides the com- 
plete answer to the problem, useful insight can be gained 
from an examination of the perturbative approach. We 
therefore consider next the case of small spin orbit cou- 
pling and will assume that the Rashba term can in- 
deed be treated as a perturbation. In order to sim- 
plify the notation, in this Section we consider a potential 
that confines the electrons in < y < W rather than 
— W/2 < y < W/2. For simplicity we will label the per- 
turbed eigenstates via the (approximate) quantum num- 
bers of the unperturbed solutions. Accordingly we write 



¥n± (k x ,y) = ip a n ± (k x ,y) + 0(a) , 



(20) 



where 



LW W 1 )y ' 1 ' 



with n a positive integer. We also find convenient to 
choose the spinors along y, in such a way that the states 
are eigenstates of the y reflection symmetry. The relevant 
matrix elements of the spin orbit interaction are: 



{ip% ± \Ro\ip° n± )=^-%ak x 5 nn '8k x k' , 



(22) 



<^> |7? |^ ± )=T 



2ha nn'(l - (-l) n+n ') 
~W (n 2 - n' 2 ) 



We can easily obtain the form of the perturbation ex- 
pansion, up to second order: 

h 2 /irn\ 2 Ti 2 k 2 
tn±(k x )=— (—;) + 



2m \ W 2 J 2m 
Thak x -^ma 2 + 0(a 3 ) . (23) 



Higher order terms are also readily obtained. 

The necessary condition for the validity of the per- 
turbative expansion is that the matrix elements (|22|) be 
smaller than the unperturbed energy differences. The 
latter are of order „,.> . This gives: 

ft , ft 

and a < . (24) 



a < 



mk x W 2 



mW 



Clearly the perturbative expansion fails for large values 
of k x , or for a channel with large W. 

The perturbation theory allows us to explicitly calcu- 
late a number of interesting quantities. In particular the 
z-component of the spin polarization density of the con- 
fined states (see Eq. (fT9|) ) is given by: 

Pn±(fca»y) = (<Pn±(k x ,y')\a z 5(y - y')\(p n ±(k x ,y')) (25) 



otm „ , . 
T-f-F n (y/W) 



amW 



k x G n (y/W) + 0(a 3 ) 



where the functions F n and G n are given by 
64^-^' inn 

7T 



irny irmy 



G n {y/W)J^Y. 



(n 2 — m 

nm . nny . irmy 
7—5 5777 sin sin , (26 

(n 2 -m 2 ) 3 W W 



where the sum extends over positive even (odd) integers 
m for odd (even) n. 

In the same way, while 'P^{k x , y) vanishes, the corre- 
sponding perturbative result for the polarization along y 
of the n± states can be expressed as: 

Vl±{k x ,y) = (<Pn±(k x ,y')\a y S(y - y')\(p n ±(k x ,y')) (27) 



1 9R 

H n (y/W) = — 



W^W^) WH n (y/W) + 0(aP), 



= ± 

where the function H n is defined as follows 

2 ™y 



(n 2 



-TTTTT Sm 

i ) W 



m nn . nny . 7m y 
sm — ^\\\ 

i 2 — m 2 ) 2 (n 2 — n' 2 )(m 2 — n' 2 ) W W 



mm'n 2 



— ' (n 2 — m 2 ) 2 (n 2 — m' 2 ) 2 

rn I x ' ' 



. irmy . irm y 
• sm sm 



w 



w 



, (28) 



where the sums extend over the positive integers. In 
particular while nl is of the same parity as n, both m 
and m' are of opposite parity. Moreover n' ^= n. 

Notice that as expected, in all these formulas the sub- 
stitution k x — > —k x and ± — > =p gives a state with the 
same energy and opposite spin polarization. 



III. EDGE SPIN POLARIZATION 

We consider here the spin accumulation at the bound- 
aries of the wire in the presence of an electric current in 
the x direction. 
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FIG. 4: Number density N v (y) (in m~ units) of individual eigenstates at the Fermi level of a quantum wire. The parameters 
are the same as in Figure [2] The wave vector values are determined by the crossings of the bands with the horizontal dashed 
line at 6 meV in Figure [2] and are given at the top (dots on the kx line). The corresponding N v {y) functions are shown in 
the panels, where the numbers label the first twelve wave vectors, counted from left to right on the top k x line. States with 
opposite wave vectors have the same density. 




FIG. 5: Spin polarization density V„(y) along the z direction (in m^ 1 units) of the same eigenstates of Figure[4] States with 
opposite wave vectors have polarization of opposite sign. Solid lines correspond to the first twelve k x values, counted from left 
to right on the top k x line, while dashed lines refer to — k x . 



Although the response of a system of electrons sub- 
ject to Rashba spin orbit and an applied electric field 
is quite complicated and requires a careful analysis, 2,3 
for simplicity sake we will model here a current carrying 
steady state by assuming a modified electron occupation 
distribution. In particular we will assume that a finite 
chemical potential difference <5/i = /i+ — /i_ is established 
between the right moving and left moving states. This 
is depicted in Figure [5] by the two short horizontal solid 
lines (whose separation from the equilibrium Fermi level 
- dashed - is not in scale). The present model is meant 
to describe ballistic transport with all the chemical po- 



tential drop occurring at the contacts at x — ±oo. 

Within this framework, the total spin polarization den- 
sity due to the current flux can be obtained in the linear 
regime from the relation: 

M{y) = ±if-Y^-^—V v {k Fv ,y) , (29) 
1 Zirn vp v 

where the sum is restricted to the occupied one- 
dimensional subbands with kp v and v Fv being the (pos- 
itive) Fermi wave vector and Fermi velocity relative to 
the occupied band as labeled by the index v. In this for- 
mula the positive sign refers to the case of hole transport. 



6 




FIG. 6: Spin polarization density Viiy) along the y direction (in m _1 units) of the same eigenstates of Figure [4] States with 
opposite wave vectors have polarization of opposite sign. Solid lines correspond to the first twelve k x values, counted from left 
to right on the top k x line, while dashed lines refer to — k x . 



As a consequence the direction of the magnetization only 
depends on the current direction and not on the type of 
carrier involved. 

The perturbative result for the z component M. z is 
readily obtained from Eqs. (J^HJ) : 

MM * T6^(^-) 2 Y,G n (y/W), (30) 

where the sum is limited to the positive integers corre- 
sponding to the occupied bands. The results of the (ex- 
act) numerical calculations of M. z (y) are shown in Figure 

m 

Interestingly the expression given in (|27[) is not suf- 
ficient to obtain the leading corresponding perturbative 
expression for A4 y (y). The reason is simply that in the 
confined geometry the first finite contribution to this 
quantity is clearly of order a 3 (see also below). This 
term can be of course readily calculated. The results for 
the exact calculation of A4 y (y) are plotted in Figure [8] 
alongside the perturbative expression. 

It is remarkable that the z component of the magneti- 
zation is larger for small values of the spin-orbit coupling 
constant while the two components become comparable 
for larger values. Furthermore, while M. z {y) changes sign 
and displays an asymmetric behavior across the wire, 
Ai y (y) remains of the same sign and is instead symmet- 
ric. 

As expected, in all cases, there is good agreement with 
the approximate analytical formula (|30[) for small values 
of a. 

Our theory can be extended to the case of non ballistic 
transport. A simple treatment of such a situation can be 
obtained by setting Sfj, — > eEvp v T v in Eq. (|29p. Here E 
is the magnitude of the driving electric field and vf v t v 
is the elastic mean free path appropriate for the states 



of the one dimensional subband u, a length scale that 
we assume much larger than the width W of the wire. 
By taking for simplicity the same scattering time for all 
subbands r„ = r the ensuing magnetization can be ob- 
tained. The results for M z {y) are plotted in Figure [9] 
We notice that the oscillations of the magnetization are 
somewhat damped as compared to the ballistic case. On 
the other hand, the overall sign of the effect is unchanged, 
the behavior being in qualitative agreement with the ex- 
perimental results of Ref. [^. We haste to state however 
that our model is rather different from the situation of a 
dirty three-dimensional sample dealt with in the experi- 
ment. 



IV. DISCUSSION 

As we have noted, while for small spin orbit the out 
of plane (z-axis) polarization M z (y) induced by an x- 
axis current is quadratic in the coupling constant, the in 
plane (y-axis) polarization M y (y) is much smaller, first 
appearing in third order. This result is only in apparent 
contradiction with the well known corresponding bulk re- 
sult in two dimensions which entails a linear dependence 
of M. y (y) on the coupling constant. Within the simplified 
model used in Section IIIII such a behavior can be read- 
ily derived by assuming the distorted momentum space 
electron distribution function of Figure 1101 The final ex- 
pression is given by 

J^„> . °f a ; 5 , (3D 

■L ir A n y2eF + mar 

where <5/i = 2Tivp8k. As expected for small a M y (y) is 
linear in the coupling constant J£ 

The difference lies with a quenching of the y-axis po- 
larization due to the confining potential. As it can be 
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FIG. 7: Spin polarization density in the z direction of a W = 
0.1 /xm wide quantum wire, for carriers with effective mass 
m = 0.3mo and g = 2. The Fermi energy is = 6 meV, 
which corresponds to an electron density of 0.67 10 12 cm -2 . 
We used Sfi — 0.1 meV which gives a current in the wire of 
23 nA. In the top panel ha = 1.5 10 -9 meVm while in the 
bottom one ha = 3.5 10 -9 meVm. The perturbative result is 
also plotted (dashed). As a comparison of the magnitude of 
the effect, the Pauli susceptibility is \p = 72.5 /ib/(/iiii 2 T). 



surmised by inspecting Eq. (f2"T)) , this effect can be readily 
seen to originate from the lowest order cancellation of the 
contributions to M y (y) stemming from states within the 
same spin-split one dimensional band. The cancellation 
only occurs when the spin-splitting is smaller than the 
energy quantization associated with the one-dimensional 
confining potential. This is of course not operational in 
the bulk since in such a situation the energy spectrum is 
continuum. 

The crossover between these two regimes approxi- 
mately occurs when the spin-splitting equals the energy 
spacing between the eigenvalues of the confining poten- 
tial along the y direction, i.e. for m ° ff ~ 1. In order to 
exemplify this phenomenon we have plotted in Figure [TT1 
the spatially averaged magnetization density 

1 fW/2 

M y = w / M v(y) A v > ( 32 ) 

W J-W/2 

as a function of the value ha of the spin orbit coupling 
constant. With the parameters chosen in Figure [TT1 the 
crossover occurs approximately for ha = 1.7 10~ 9 meV m. 
The dashed curve superimposed to the tp = 2.5 meV 
dotted line is purely cubic and provides a useful guide to 
the eye for that case. 



FIG. 8: Spin polarization in the y direction corresponding 
to the same parameters of the previous Figure [7] In the 
top panel ha — 1.5 10~ 9 meVm while in the bottom one 
ha = 3.5 10 -9 meVm. The leading perturbative contribu- 
tion (cubic in a) is also plotted (dashed). 
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FIG. 9: Magnetization density of a W — 0.1 /J,m wide quan- 
tum wire, for carriers with same parameters as the first panel 
of Figure [7] We used E — 1 kV/m and a constant value for 
r = 7.3 10 -13 s, which gives the same current in the wire of 
Figure [7] The perturbative result is also plotted (dashed) . 
The magnetization is somewhat different from the previous 
case, where the oscillations of the magnetization are more 
pronounced. The overall sign of the effect is the same. 



Exactly the opposite mechanism is responsible for the 
case of M z (i/), which, as we have seen, is instead vanish- 
ing in the bulk. Interestingly, as one can infer from Figure 
[TT1 the expression of Eq. (|3ip . formally valid in the bulk 
for finite values of a when in principle the perturbative 
condition is inapplicable, still gives good results for the 
case of a rather narrow channel with W = 0.15 /zm. 
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FIG. 10: Distorted momentum space occupation distribution 
used to calculate the bulk value of the magnetization. K± are 
defined in Eq. ([5]). 




FIG. 11: Average spin polarization density in the y direction 
as a function of a. The parameters used are: W = 0.15 /im, 
€f = 2.5, 10, 40 meV, <5/i = 0.05 meV, and m — 0.3mo (mass 
as in previous Figures). Here the dots represent the (exact) 
numerical result, while the dotted curves are obtained from 
Eq. (|3ip in the text. The dashed line is a guide to the eye 
representing a pure cubic behavior. 



cussed are compatible with a quite general and rigorous 
argument regarding the parity of M. z (y) and M y {y) with 
respect to the sign of the spin orbit coupling constant. 
Consider the effect on the hamiltonian (including contri- 
butions from the electron electron interaction and non 
magnetic impurities) of the application of the operator 
a z . Clearly a z H(a)a z — H(—a). Then it immediately 
follows that given an eigenstate of the problem for a given 
value of a, the corresponding solution for —a can be ob- 
tained simply by application of a z . This in turn implies 
that, quite generally, while M. z (y) is even with respect 
to a, M. y {y) is odd. 

We have shown that a net spin polarization is estab- 
lished when a steady state current is run through the wire 
obtained by laterally confining a two dimensional elec- 
tron gas in the presence of Rashba spin orbit. Reversing 
the current inverts the spin polarization. The effect is 
phenomenologically analogous to, yet distinct from the 
spin-Hall effect that has recently been the subject of in- 
tense discussion (for a useful review, see Ref. [Ill) and 
was demonstrated experimentallyi£ i 12 ' 13 i 14 The mecha- 
nism relevant to the present discussion simply stems from 
the structure of the electronic wavefunctions in a con- 
fined geometry. Accordingly, there is no need to make 
appeal to spin currents or impurity scattering. More- 
over its geometric nature should leave the phenomenon 
mostly unchanged in the presence of a moderate amount 
of impurities in the wire. A similar effect should also be 
relevant in the case of electronic states localized by an 
impurity potential as well as in the transition to a low 
density Wigner crystal in the presence of Rashba spin 
orbit. 
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